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Abstract. The electron transport through the triple-barrier resonant tunnelling 
diode (TBRTD) have been studied by the self-consistent numerical method for the 
Wigner-Poisson problem. The electron flow through the TBRTD can be controlled by 
the gate voltage applied to one of the potential well regions. For different gate voltage 
values we have determined the current-voltage characteristics, potential energy profiles, 
and electron density distribution. We have found the enhancement of the peak-to- valley 
ratio (up to ~10), the appearance of the linear current versus bias voltage behaviour 
within the negative-differential resistance region, and the bistability of the current- 
voltage characteristics. The analysis of the self-consistent potential energy profiles and 
electron density distribution allowed us to provide a physical interpretation of these 
properties. 



1. Introduction 

The standard analysis of transport processes in solids is based on the semiclassical 
Boltzmann transport equation, e.g. [H, El EJ H]. However, in the solid-state 
nanostructures, the quantum effects become more pronounced and the semiclassical 
description fails [5|. In particular, the theoretical description of electron transport in 
semiconductor nanostructures and nanodevices requires the fully quantum-mechanical 
treatment. In the literature, one can find several quantum-mechanical approaches to the 
electronic transport in nanodevices, such as those based on the envelope function [SI BE], 
density matrix 0, [101 El EG] , non-equilibrium Green function [131 HS1 IS] , and non- 
classical distribution function [I7J [TEl [19]. One of the most transparent and powerful 
methods, used to investigate the quantum electron dynamics in the nanodevices, is based 
on the quantum kinetic equation for the Wigner distribution function (WDF) [201 I2T1 122] . 

The quantum kinetic equation for the WDF has been applied to a description 
of the electron transport in the resonant tunnelling diode (RTD), i.e. the double 
barrier semiconductor nanostructure [UJ [231 [Ml [251 I2S]- The triple-barrier resonant 
tunnelling diode (TBRTD) [271 I2H] is a natural extension of double barrier RTD. In 
the TBRTD, we deal with two quantum wells separated by the barrier regions. The 
resonant electron tunnelling, which is responsible for the successful operation of the 
resonant-tunnelling nanodevices, strongly depends on the relative positions of energy 
levels of the electrons confined in the different quantum-well regions. Therefore, the 
TBRTD nanodevice provides a possibility of additional tuning of the resonant condition 
by applying a gate voltage to the one of the quantum-well regions. The preliminary 
results for the TBRTD nanodevice are presented in our previous paper [29J. These 
results |29j have been obtained using the model that neglects the electrostatic field 
produced by the electrons gathered in the nanodevice. In the present paper, we study 
the electron transport through the gate-voltage controlled TBRTD nanodevice taking 
into account this electrostatic field by a self-consistent procedure. For this purpose, we 
apply a numerical method of self-consistent solution of the quantum kinetic equation 
for WDF and the Poisson equation. As a main result, we obtain the current-voltage 
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characteristics of the TBRTD as functions of nanodevice parameters and gate potential. 
This approach allows us to simulate the operation of the TBRTD nanodevice. 

The present paper is organized as follows. In Sec. [2] we introduce the quantum 
kinetic approach and describe the model of the nanodevice. In Sec. [3] we present the 
results and their discussions, and in Sec. HI we provide the conclusions and summary. 

2. Theoretical Method and Model of Nanodevice 

The resonant-tunnelling devices consist of different semiconductor layers which act as 
potential barriers and wells for the charge carriers. The transport properties of RTD and 
TBRTD are determined by the chemical and physical properties of the semiconductors 
forming the layers as well as by the thicknesses of the layers. At the semiconductor 
heterojunction, the difference in energy gaps between two different materials leads to 
a conduction (valence) band discontinuity. The position-dependent conduction band 
minimum forms the potential profile for electrons with potential barrier and well 
regions. In the potential well regions, the electrons can form the quasi-bound states. 
These states strongly affect the electron transport via the nanostructure due to the 
quantum interference effects which - under certain conditions - can lead to the resonant 
tunnelling. If the layers are fabricated from the doped semiconductors, the thermal and 
electric-field induced ionization of donor and acceptor impurities yields additional charge 
carriers, which are non-uniformly distributed over the nanostructure. These charge 
carriers are sources of an additional electric field which leads to the second component 
of the electron potential. The theoretical description should take into account both the 
components of the electron potential energy, whereas the density of the additional charge 
carriers and its electrostatic potential should be calculated by a self-consistent manner. 
The present approach is based on the self-consistent solution of the quantum kinetic 
equation for the Wigner distribution function and Poisson equation for the electrostatic 
potential. 

In the quantum kinetic theory, the conduction band electrons are described by the 
WDF. If we take on the x axis in the growth direction of the semiconductor layers 
forming the TBRTD and assume the translational symmetry in the lateral directions 
(y,z) (cf. figured]), we can obtain the basic electron transport characteristics of the 
TBRTD from the Wigner distribution function p w (x,k,t). Neglecting the scattering 
processes and intervalley transitions in the conduction band, we write the quantum 
kinetic equation in the form [30, ETJ 

op w (x,k,t) %k op w [x,k,t) i f , , , , , , . 

ST 1 ^ = / dk U w( x i k ~ k )Pw{X, k , t) , (1) 

at m ox 2nn J 

where m is the conduction band effective mass and the integral kernel U w (x,k — k') 
represents the non-local potential energy defined as 

U w (x, k-k') = J dx' [U(x + x'/2) - U(x - x'/2)} exp [ - i(k - k')x'} , (2) 

where U(x) is the total electrostatic potential energy of the electron. 



Self- consistent Wigner distribution function study of gate-voltage controlled TBRTD 4 



In present work, we consider the steady-state solutions of ([Tj), i.e., the solutions 
obtained under assumption dp w (x,k,t)/dt = 0. The potential energy U(x) consists 
of the conduction band potential Ub{x) and the Hartree potential energy Uh{x) of 
the electrons which originates from the ionized donor impurities (the n-type doping is 
assumed). Explicitly, 

U(x) = U B (x) + U H (x). (3) 

The conduction band component Ub{x) can be written down in the form 

N 

U B (x) =Y,Ui<S)(x - Xi)Q(x i+1 - x) , (4) 
i=i 

where N is the number of barrier regions, Xi is the position barrier-well (well-barrier) 
interface, i.e., — x,| determines the thickness of the barrier region, Q(x — Xj) is the 
Heaviside step function, and [/, is the height of the z-th barrier. 

The Poisson equation for the Hartree potential energy takes on the form 

d 2 U H (x) e 2 r , T , . . 

, 2 =—[N D (x)-n(x)] , (5) 
dx z e e 

where e is the elementary charge, Eq is the vacuum electric permittivity, e is the relative 
static electric permittivity, Nd(x) is the concentration of the ionized donors, and n(x) 
is the electron density. 

If we determine the WDF, we can find from its zero and first moments the electron 
density and current density, respectively [30]. In the explicit form, 

n ( x ) = 7T f dkp w (x,k) (6) 

and 

j( x ) = TT [ dk—p w (x,k). (7) 

The quantum kinetic equation (JTJ and Poisson equation (jSJ) should be solved 
simultaneously by a self-consistent procedure [32]. We assume the Dirichlet boundary 
conditions for the Poisson equation: Uh(0) = and Uh(L) = —eVb, where the bias 
voltage Vb = Vr — Vl is applied between the right (R) and left (L) electrodes separated 
by distance L. Moreover, for the WDF we take on the open boundary conditions which 
have the form [33] 



Pw(0, k) 
Pw(L,k) 



= f L (k) , (8) 

fc>0 

= f R {k) , 

fc<0 

where f L< ^ R \k) is the supply function [31], i.e., the Fermi-Dirac distribution function 
integrated over the transverse momenta of electrons in the left (right) reservoirs. The 
supply function for the left (right) reservoir has the form 

2 7„2 
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Figure 1. Self-consistent potential energy profile in the TBRTD for zero bias voltage. 
The position is measured in the layer growth direction, (J.l(r) is the electrochemical 
potential of the left (right) contact, Uo is the height of the potential barrier, and U g is 
the gate- voltage controlled shift of the bottom of the right potential well. The undoped 
region consists of three Alo.3Gao.7As potential barriers and two GaAs potential wells 
sandwiched between two spacer regions. 



where T is the temperature and i^lir) is the electrochemical potential of the left (right) 
reservoir. 

In the present paper, we consider the TBRTD nanodevice composed of GaAs 
and Alo.3Gao.7As layers (figure CD). The TBRTD nanodevice consists of the triple 
barrier region with two potential wells, spacer layers, and contacts (cf. figure [I]). 
The contacts are made from the n-doped GaAs with homogeneous donor concentration 
No = 2 x 10 18 cm -3 . The calculations have been performed for the following parameters 
of the layers: the thickness of each contact is equal to 17 nm, the thickness of each spacer 
layer is 3 nm, the thickness of each potential barrier (well) is 3 nm (5 nm), the height 
of the potential barrier Uo = 0.27 eV, the total length of the nanodevice is L = 59 nm. 
In the calculations, we change energy U g of the bottom of the right potential well which 
allows us to determine the effect of gate voltage applied to this region. We assume 
that the electrons are described by the conduction band effective mass of GaAs, i.e., 
m = 0.667 mo, where mo is the free electron mass. Similarly, we take on the electric 
permittivity, e = 12.9, and the lattice constant, a = 0.565 nm for GaAs. 

Temperature T enters into the problem via supply function ([9]), which exhibits an 
almost unchanged shape in the interval < T < 77 K. However, taking on the larger 
value of T considerably accelerates the calculations. Therefore, in the present paper, we 
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simulate the low-temperature properties of the TBRTD putting T = 77 K. 

We have obtained the steady-state solutions of the Wigner- Poisson problem for the 
TBRTD using the numerical scheme proposed by Biegel and Plummer [M] and Kim |35j . 
In the first step, we find the steady-state Wigner function using the time-independent 
form of (JTJ) with the potential energy given by (Tj0) and calculate the electron density 
according to ([6]). In the next step, after inserting the electron density into (jSJ) we 
calculate the new potential energy profile. Using this potential energy, we again solve 
the time-independent form of (JTJ). This procedure is repeated until the convergence 
is reached. All the simulations were carried out using the computational grid with 
N x = 106 mesh points for the position x and Nk = 72 for the wave vector k. 

3. Results 

In the gate-controlled TBRTD nanodevice, we have a possibility of designed changing 
the energy levels of the electrons confined in the quantum wells. In the considered 
TBRTD model, we can shift up and down the energy levels in the right potential well 
by changing energy U g of the potential well bottom. In a practical realization, U g 
can be changed by applying gate voltage V g to the nearby gate attached to the right 
potential well layer. Then, U g = —eaV g , where a is the voltage-to-energy conversion 
factor which is characteristic to the nanodevice [36] . Alternatively, U g can be changed in 
the technological process by fabricating the right-potential well layer from Al^Gai-^As 
with x < 0.3. 

The resonant tunnelling devices are characterized by the Peak-to- Valley Ratio 
(PVR) which is defined as the ratio of the maximum (peak) current to the minimum 
(valley) current on the current-voltage characteristics. The results of our previous 
calculations [29] suggest the large value of PVR (PVR = 27) obtained for U g = 0.1 eV. 
The present more realistic approach leads to smaller values of PVR (figure [2]). The 
smaller values of the PVR result from the considerable modification of the potential 
profile due to the electrostatic interaction with the localized electrons (cf. figure [p. 
The maximum value of the PVR obtained in the present calculations is 9.9. We note, 
however, that this value of the PVR is still larger than those determined for the resonant- 
tunnelling diodes based on the GaAs technology [37]. 

Figure [2] shows the current-voltage characteristics of the TBRTD nanodevice 
calculated for different potential energy bottom U g of the right potential well. We 
observe that the increasing U g shifts the positions of the current peaks towards the 
higher bias voltage. This effect results from the increasing ground-state energy of the 
electron in the double well potential. Then, the condition of resonant tunnelling is 
satisfied at the higher bias voltage. We have found that for U g < 0.13 eV the PVR 
increases with the increasing U g . At U g = 0.13 eV the PVR reaches the value 9.9 which 
is the maximal value of the PVR obtained in the present calculations. 

The obtained features of current-voltage characteristics can be explained if we 
consider the self-consistent potential energy profiles and electron distribution in the 
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Figure 2. Current -voltage characteristics and PVR coefficients of the TBRTD 
nanodevice calculated for different U g . 

TBRTD. Figure [3] displays the results for U g = 0.09 eV. We observe that the electron 
distribution for the bias corresponding to the maximum of current (VI = 0.31 V) 
considerably differs from those for the minimum of current (H = 0.44 V). The resonant 
tunnelling (current maximum) appears for the large electron accumulation in the left 
potential well which results from the fact that - in this case - the tunnelling current 
flows via the one-electron ground state which is localized in this potential well. For the 
minimal current the electrons occupy both quantum wells with the comparable density. 

A closer inspection of the current- voltage characteristics for U g = 0.13 eV (figure [2]) 
exhibits a notch in the regime of negative differential resistance. If the bias voltage 
becomes larger than 0.41 V, the current density stops to decrease but starts to increase 
linearly with the increasing bias. This increase occurs in the narrow bias voltage regime 
of the width estimated to be 0.02 V. This linear increase of the current is sometimes 
called the "plateau-like" behaviour [3E]. For VI exceeding 0.43 V we again observe the 
negative differential resistance (cf. figure [2]). This unusual property of the current- 
voltage characteristics results from the formation of the potential well in the region of 
the left spacer, which appears in the narrow bias voltage regime (figure H]). The electron 
density distribution and potential profile in the bias voltage regime 0.41 V < VI < 
0.43 V are displayed in figure HI In this spacer-related potential well, the electrons can 
form the quasi-bound states with the discrete energy levels. The ground-state energy 
of the electron in the left spacer is slightly lower than the energies of the quasi-bound 
states in the quantum wells which are responsible for the resonant tunnelling through 
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Figure 3. Potential energy profile and electron density distribution for U g = 0.09 eV 
and bias voltage V& = 0.3f V and 0.44 V that corresponds to the current maximum 
and minimum, respectively. 

the TBRTD. The spacer-related potential well leads to the accumulation of electrons 
flowing from the left contact. Because the difference between energy levels in the left 
contact and in the left quantum well is very small, the conditions for resonant tunnelling 
are approximately fulfilled. In the bias voltage regime 0.41V^ < Vj, < 0A3V the electron 
can flow from the left contact via the spacer and the triple barrier region to the right 
contact. Therefore, the notch on the current-voltage characteristics which appears in 
this bias voltage regime, is a signature of the resonant tunnelling through the spacer 
quasi-bound states. If the bias voltage increases above Vb = 0.43 V, the spacer-related 
potential well becomes more shallow and finally disappears which rises the energy of the 
electrons in the spacer. This leads to the charge outflow from the left spacer (cf. the 
curves for VJ, = 0.43 V in figure H]). At higher bias voltage the resonant tunnelling is 
no longer responsible for the electron transport which leads to an abrupt decline of the 
current. 

The change of the current- volt age characteristics which appears for U g > 0.13 eV, 
results from the considerable modification of the potential profile which is visible if we 
compare the potential profile for U g = 0.13 eV (figure H]) with those for U g = 0.09 eV 
(figure [3]) . The results presented in figure [3] show that the spacer-related potential well 
almost disappears. Therefore, there is no charge accumulation in the spacer region. The 
plots presented in figure [3] correspond to the usually observed behaviour of current vs 
voltage in the RTD. 

If U g exceeds 0.13 eV, the negative-differential resistance behaviour is further 
modified (figure [6]): the notch becomes more distinct and the width of the linear- 
increase region is extended. At U g = 0.15 eV, the notch dominates in the current- 
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Figure 4. Potential energy profile and electron density distribution for U g — 0.f3 eV 
and bias voltage V b = 0.41, 0.42, 0.43 V. 

voltage characteristics and the current maximum has the form of the sharp peak reached 
after the linear-increase behaviour [figure [01(c)] . The linear current versus bias voltage 
behaviour has been obtained for the increasing as well decreasing bias voltage. For the 
decreasing bias voltage and U g = 0.15 eV we have obtained the following two linear 
current-versus-voltage regimes: (i) 0.32 V < V& < 0.36 V, and (ii) 0.37 V < < 
0.38 V [cf. bold lines inline)]. The potential profile and electron density for these two 
regimes are depicted in figure El The difference in the electron distribution (figure [5]) for 
these two regimes results from the formation of quasi-bound states in the spacer layer. 
This in turn leads to the appearance of distinct regimes (i) and (ii). In regime (i), the 
transmission via the spacer-related ground state dominates in the resonant tunnelling, 
while in regime (ii), we deal with the tunnelling via the first excited state. We have 
checked the above conclusions by a direct numerical calculation of the quasi-bound state 
energy levels. 

Additionally, we have found that the current-voltage characteristics exhibit the 
bistability (figure [6]). Figures [6|a, b, c) show that the current- volt age characteristics 
obtained for the increasing bias (solid curves) differ from those for the decreasing bias. 
The largest difference is observed in the region of maximal current. 

The bistability in TBRTD current- voltage characteristics can be explained by either 
the mathematical or physical arguments. If we consider the mathematical properties of 
the self-consistent Wigner-Poisson problem, we observe its nonlinearity. The solutions 
of the nonlinear problem are very sensitive to the initial conditions, therefore, they 
depend on the direction of the bias voltage change. Physical there are two orgins of 
bistability. First of them is the feedback between the charge accumulated in quantum 
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Figure 5. Potential energy profile and electron density distribution for U g — O.f 5 eV 
and for bias voltages V& = 0.32 V and 0.38 V, which correspond to the two linear-growth 
regimes in the current- voltage characteristics depicted in figure [6] (c). 



wells and the flowing current [391 SO], and the second is a current oscillation [41] . The 
bias voltage V h = 0.41 V corresponds to one of the two bistability points in the current- 
voltage characteristics as it is shown in figure [3 We observe that both the potential 
energy and electron density are considerably different in the cases of increasing and 
decreasing bias voltage. The increasing bias causes that the electrons are localized in 
both the potential wells, however, the electron localization in the left quantum well 
is stronger (cf. solid curves in figure [7]). The electron localization results from the 
formation of the quasi-bound states in both the potential well regions, which - at the 
resonant conditions - leads to the increase of the current. This type of behaviour is 
typical for asymmetric potential profile. On the contrary, the decreasing bias leads to 
the accumulation of electrons in the left spacer (near the left potential barrier). In this 
case, the electron density is almost zero within the potential wells which results in the 
considerable reduction of the maximal current [cf. figure EJ^c)]. 

4. Conclusions and Summary 

In the present paper, we have studied the transport properties of electrons in the triple- 
barrier resonant tunnelling diode. We have solved the Wigner-Poisson problem by the 
self-consistent procedure and determined the current- voltage characteristics for different 
gate voltage applied to the right potential well region. We have found the fairly large 
value (~10) of the peak-to- valley ratio. The enlargement of the peak-to- valley ratio over 
that observed for the double-barrier resonant tunnelling diodes is caused by appearance 
of the linear regime within the negative differential resistance region that is produced 
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Figure 6. Current-voltage characteristics of the TBRTD for (a) U g = 0.13 cV, (b) 
U g = 0.14 cV, and (c) U g = 0.15 eV. Solid (dashed) curves with arrows down (up) 
correspond to the results obtained for increasing (decreasing) bias voltage. 



by the join effect of the additional potential barrier and potential well. The additional 
potential barrier leads to the accumulation of electrons in the left spacer, while the 
second potential well leads to the formation of quasi-bound states in this potential well. 
For properly chosen nanodevice parameters the resonance conditions are satisfied for 
the quasi-bound states appearing in the three potential well regions, namely, the left 
spacer and two potential wells. This leads to a considerable increase of the resonant 
current peak. 

We have also found interesting properties of the current- volt age characteristics of 
the TBRTD: the change of the negative differential resistance into the positive one (the 
occurrence of the notch in the regime of the linear dependence of the current on the 
bias voltage) and the bistability. We have explained these properties by analyzing the 
corresponding changes of the potential profile and electron density which are controlled 
by the gate potential. In the present study, we focus on the influence of the self- 
consistent Hartree potential on the transport properties of the TBRTD, therefore, we 
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Figure 7. Potential energy profile and electron density distribution in the TBRTD for 
U g = 0.15 eV and for bias 14 = 0.41 V which corresponds to the bistability region in 
figureEl Solid (dashed) curves show the results obtained for the increasing (decreasing) 
bias voltage. 



neglect the electron-phonon scattering. The bistable behaviour found in the current- 
voltage characteristics of the TBRTD requires a further study that should include the 
time dependence of the current. 

In summary, the current-voltage characteristics of the triple-barrier resonant 
tunnelling diode shows a variety of interesting properties that can be effectively tuned 
by applying the gate voltage. 
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